Weak Galerkin Finite Element Method Research Articles

Weak Galerkin finite element methods for semilinear Klein–Gordon equation on polygonal meshes

Weak Galerkin finite element methods for semilinear Klein–Gordon equation on polygonal meshes

Optimal [formula omitted] error estimates of stabilizer-free weak Galerkin finite element method for the drift-diffusion problem

We continue our effort in Li et al. (2024) to develop a stabilizer-free weak Galerkin finite element method (SFWG-FEM) for the two-dimensional unsteady-state drift-diffusion (DD) model. The transient state equation includes both the first derivative convection and the second derivative diffusion terms, and it couples with a Poisson potential equation. The piecewise Pk(k≥1) polynomials elements in the interior and on the boundaries are utilized for both electron concentration and electric potential functions. The vector function of [Pj(K)]2(j≥j0) for the discrete weak gradient space is used on each element K. The main difficulty is the treatment of the non-linearity and coupling of the system. Based on an introduced Ritz projection and two L2 projection operators, optimal error estimates of O(Δt+hk+1) in the L2 norm and O(Δt+hk) in the energy-like norm are derived. Numerical experiments are presented to illustrate the theoretical analysis.

A Stabilizer Free Weak Galerkin Finite Element Method for Elliptic Equation with Lower Regularity

A Stabilizer Free Weak Galerkin Finite Element Method for Elliptic Equation with Lower Regularity

Analysis and Computation of a Weak Galerkin Scheme for Solving the 2D/3D Stationary Stokes Interface Problems with High-Order Elements

Abstract In this paper, we present a high-order weak Galerkin finite element method (WG-FEM) for solving the stationary Stokes interface problems with discontinuous velocity and pressure in ℝ d (d = 2, 3). This WG method is equipped with stable finite elements consisting of usual polynomials of degree k ≥ 1 for the velocity and polynomials of degree k – 1 for the pressure, both are discontinuous. Optimal convergence rates of order k + 1 for the velocity and order k for the pressure are established in L 2-norm on hybrid meshes. Numerical experiments verify the expected order of accuracy for both two-dimensional and three-dimensional examples. Moreover, numerically it is shown that the proposed WG algorithm is able to accommodate geometrically complicated and very irregular interfaces having sharp edges, cusps, and tips.

A modified weak Galerkin finite element method for the Maxwell equations on polyhedral meshes

We introduce a new numerical method for solving time-harmonic Maxwell’s equations via the modified weak Galerkin technique. The inter-element functions of the weak Galerkin finite elements are replaced by the average of the two discontinuous polynomial functions on the two sides of the polygon, in the modified weak Galerkin (MWG) finite element method. With the dependent inter-element functions, the weak curl and the weak gradient are defined directly on totally discontinuous polynomials. Optimal-order convergence of the method is proved. Numerical examples confirm the theory and show effectiveness of the modified weak Galerkin method over the existing methods.

A posteriori error estimate of the weak Galerkin finite element method solving the Stokes problems on polytopal meshes

AbstractIn this paper, we propose an a posteriori error estimate of the weak Galerkin finite element method (WG‐FEM) solving the Stokes problems with variable coefficients. Its error estimator, based on the property of Stokes' law conservation, Helmholtz decomposition and bubble functions, yields global upper bound and local lower bound for the approximation error of the WG‐FEM. Error analysis is proved to be valid under the mesh assumptions of the WG‐FEM and the way can be extended to other FEMs with the property of Stokes' law conservation, for example, discontinuous Galerkin (DG) FEMs. Finally, we verify the performance of error estimator by performing a few numerical examples.

A Crank-Nicolson WG-FEM for unsteady 2D convection-diffusion equation with nonlinear reaction term on layer adapted mesh

This paper is contributed to explore how a Crank-Nicolson weak Galerkin finite element method (WG-FEM) addresses the singularly perturbed unsteady convection-diffusion equation with a nonlinear reaction term in 2D. The problem and some asymptotic behavior results are given for the exact solution and its derivatives with the parameter ε. These results are essential for proving the uniform convergence of the proposed WG-FEM. A tensor product Shishkin mesh featuring piece-wise discontinuous bilinear polynomials is employed to handle the layers for uniform convergence at different parameter values of ε. An error estimate ‖uN−INu‖WG is presented, where INu denotes the vertex-edge-cell interpolation of the solution u, and ‖⋅‖WG denotes the Weak Galerkin norm. The optimal uniform convergence order is demonstrated for semi-discrete and fully discrete schemes. Various numerical experiments are conducted to validate the optimal order of convergence demonstrated by the proposed method.

The lowest-order weak Galerkin finite element method for linear elasticity problems on convex polygonal grids

This paper presents the lowest-order weak Galerkin finite element method for linear elasticity problems on the convex polygonal meshes. This method uses piecewise constant vector-valued spaces on element interiors and edges. The discrete weak gradient space introduced by this paper is the matrix version of CW0 space. The discrete weak divergence space is piecewise constant space on each element. This method is simple, efficient, stabilizer-free and symmetric positive-definite. The optimal error estimates in discrete H1 and L2 norms are presented. Numerical results are given to demonstrate the efficiency of algorithm and the locking-free property.

A robust weak Galerkin finite element method for two parameter singularly perturbed parabolic problems on nonuniform meshes

In this study, we proposed a weak Galerkin finite element method (WG-FEM) for solving two-parameter singularly perturbed parabolic problems (TP-SPPPs) of convection–diffusion–reaction type on nonuniform mesh. The WG-FEM approach incorporates an interpolation operator IN and achieves optimal order convergence in an energy-like norm for continuous and fully discrete schemes. In the fully-discrete analysis, we combined WG-FEM for spatial space on a Shishkin-type mesh with the Crank–Nicolson scheme for temporal discretization on an equidistant mesh. Since the energy norm is insufficient to capture the behavior of the boundary layer functions accurately, we have derived the optimal order of convergence in a strong balanced norm. This norm exhibits an order of O(N−1lnN)k in spatial convergence and second-order convergence in time when using an L2-projection QN as an intermediate operator. We conducted numerical examples to validate the proposed method’s theoretical findings. The results of these experiments confirmed the theoretical conclusions and demonstrated the robustness of the proposed method.

A posteriori error estimates of the weak Galerkin finite element methods for parabolic problems

In this paper, We propose the residual-based a posteriori error estimator of the weak Galerkin finite element method with the backward Euler time discretization for the linear parabolic partial differential equation. For the a posteriori error estimator, we introduce the Helmholtz decomposition technique to prove its reliability. We mainly study WG element (Pj(K),Pℓ(∂K),V(K,r)=RTj(K)). Numerical experiments based on the lowest order case, i.e. (P0(K),P0(∂K),RT0(K)), are provided to verify the theoretical research.

Generalized weak Galerkin finite element methods for second order elliptic problems

This article proposes and analyzes the generalized weak Galerkin (gWG) finite element method for the second order elliptic problem. A generalized discrete weak gradient operator is introduced in the weak Galerkin framework so that the gWG methods would not only allow arbitrary combinations of piecewise polynomials defined in the interior and on the boundary of each local finite element, but also work on general polytopal partitions. Error estimates are established for the corresponding numerical functions in the energy norm and the usual L2 norm. A series of numerical experiments are presented to demonstrate the performance of the newly proposed gWG method.

The uniform convergence of a weak Galerkin finite element method in the balanced norm for reaction–diffusion equation

In this paper, a weak Galerkin finite element method is implemented to solve the one-dimensional singularly perturbed reaction–diffusion equation. This weak Galerkin finite element scheme uses piecewise polynomial with degree k≥1 in the interior part of each element and piecewise constant function at the nodes of each element. The existence and uniqueness of the weak Galerkin finite element solution are proved. Based on the interpolation operator and the corresponding approximation properties, an ɛ−uniform error bound of O(N−(k+1)+ɛ(N−1lnN)k) in the energy-like norm is investigated rigorously. Furthermore, an ɛ−uniform error bound of O(N−(k+1/2)+(N−1lnN)k) in the balanced norm is established by the weighted local L2 projection and its corresponding approximation properties. Finally, numerical experiments validate the theoretical results. Moreover, numerical results show that this weak Galerkin finite element solution achieves the convergence rate of O((N−1lnN)k+1) in the L2 norm and the discrete L∞ norm uniformly with respect to the singular perturbation parameter ɛ.

A discrete-ordinate weak Galerkin method for radiative transfer equation

This research article discusses a numerical solution of the radiative transfer equation based on the weak Galerkin finite element method. We discretize the angular variable by means of the discrete-ordinate method. Then the resulting semi-discrete hyperbolic system is approximated using the weak Galerkin method. The stability result for the proposed numerical method is devised. A priori error analysis is established under the suitable norm. In order to examine the theoretical results, numerical experiments are carried out.

Optimal convergence analysis of weak Galerkin finite element methods for parabolic equations with lower regularity

Optimal convergence analysis of weak Galerkin finite element methods for parabolic equations with lower regularity

High order Morley elements for biharmonic equations on polytopal partitions

This paper introduces an extension of the Morley element for approximating solutions to biharmonic equations. Traditionally limited to piecewise quadratic polynomials on triangular elements, the extension leverages weak Galerkin finite element methods to accommodate higher degrees of polynomials and the flexibility of general polytopal elements. By utilizing the Schur complement of the weak Galerkin method, the extension allows for fewest local degrees of freedom while maintaining sufficient accuracy and stability for the numerical solutions. The numerical scheme incorporates locally constructed weak tangential derivatives and weak second order partial derivatives, resulting in an accurate approximation of the biharmonic equation. Optimal order error estimates in both a discrete H2 norm and the usual L2 norm are established to assess the accuracy of the numerical approximation. Additionally, numerical results are presented to validate the developed theory and demonstrate the effectiveness of the proposed extension.

Robust globally divergence-free Weak Galerkin finite element method for incompressible Magnetohydrodynamics flow

This paper develops a weak Galerkin (WG) finite element method of arbitrary order for the steady incompressible Magnetohydrodynamics equations. The WG scheme uses piecewise polynomials of degrees k(k≥1),k,k−1 and k−1 respectively for the approximations of the velocity, the magnetic field, the pressure, and the magnetic pseudo-pressure in the interior of elements, and uses piecewise polynomials of degree k for their numerical traces on the interfaces of elements. The method is shown to yield globally divergence-free approximations of the velocity and magnetic fields. We give existence and uniqueness results for the discrete scheme and derive optimal a priori error estimates. We also present a convergent linearized iterative algorithm. Numerical experiments are provided to verify the obtained theoretical results.

Preconditioned weak Galerkin finite element method for Poisson equation by least squares reconstruction

Preconditioned weak Galerkin finite element method for Poisson equation by least squares reconstruction

Comparative Study of Weak Galerkin and Discontinuous Galerkin Finite Element Methods

Comparative Study of Weak Galerkin and Discontinuous Galerkin Finite Element Methods

Weak Galerkin method for the Navier-Stokes equation with nonlinear damping term

<abstract><p>The primary focus of this research was to investigate the weak Galerkin (WG) finite element method for the Navier-Stokes equations with damping. We established the weak Galerkin finite element numerical scheme and demonstrated the existence and uniqueness of the weak Galerkin numerical solution. Additionally, optimal errors estimates for the velocity and pressure were obtained. Eventually, numerous numerical examples were reported to validate the theoretical analysis.</p></abstract>

A weak Galerkin finite element method for singularly perturbed problems with two small parameters on Bakhvalov-type meshes

A weak Galerkin finite element method for singularly perturbed problems with two small parameters on Bakhvalov-type meshes